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GLOBAL EXISTENCE AND INCOMPRESSIBLE LIMIT IN CRITICAL 
SPACES FOR COMPRESSIBLE FLOW OF LIQUID CRYSTALS 


QUNYI BIE, HAIBO CUI, QIRU WANG* *, AND ZHENG-AN YAO 

Abstract. The Cauchy problem for the compressible flow of nematic liquid crystals in the 
framework of critical spaces is considered. We first establish the existence and uniqueness of 
global solutions provided that the initial data are close to some equilibrium states. This result 
improves the work by Hu and Wu [SIAM J. Math. Anal., 45 (2013), pp. 2678-2699] through 
relaxing the regularity requirement of the initial data in terms of the director field. We then 
consider the incompressible limit problem for ill prepared initial data. We prove that as the 
Mach number tends to zero, the global solution to the compressible flow of liquid crystals 
converges to the solution to the corresponding incompressible model in some function spaces. 
Moreover, the accurate converge rates are obtained. 
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1. Introduction and main results 

In this paper we consider the global well-posedness and incompressible limit to the following 
compressible flow of nematic liquid crystals in critical spaces: 

d t p + div(pu) = 0, 

dt(pu) + div(pu © u) — /iAu — (p + A)Vdivu + VP(p) 

< , x (1.1) 

= -£div(Vd© Vd- l|Vd| 2 Ij, 

k d t d + u • Vd = #(Ad + |Vd| 2 d), 

where p £ R is the density function of the fluid, u £ M. n (N > 2) is the velocity, and d £ 
represents the director field for the averaged macroscopic molecular orientations. The scalar 
function P € M is the pressure, which is an increasing and convex function in p. We denote 
by A and p the two Lame coefficients of the fluid, which are constant and satisfy p > 0 and 
v := A + 2p > 0. Such a condition ensures ellipticity for the operator pA + (A + p)V div and is 
satisfied in the physical cases. The constants £ > 0, 6 > 0 stand for the competition between 
the kinetic energy and the potential energy, and the microscopic elastic relaxation time (or the 
Debroah number) for the molecular orientation field, respectively. The symbol © denotes the 
Kronecker tensor product such that u © u = (ujUj)i<ij<jv and the term Vd © Vd denotes a 
matrix whose th entry is d Xi d ■ d Xj d (1 < i,j < N). I is the N x N identity matrix. To 

complete the system (tm the initial data are given by 

p\t=o = Po{x), u| t=0 = u 0 (x), d| t=0 = d 0 (a), with d 0 € S^ -1 . (1.2) 

The hydrodynamic theory of liquid crystals was first proposed by Ericksen mm2! and Leslie 
pH] in 1960s. In 1989, Lin [ 25] first derived a simplified Ericksen-Leslie equation modeling liquid 
crystal flows when the fluid is an incompressible and viscous fluid. Subsequently, Lin and Liu 
m showed the global existence of weak solutions and smooth solutions for the approximation 
system. Recently, Hong m and Lin et al. [26] showed independently the global existence of a 
weak solution of an incompressible model of system (11.R in two-dimensional space. Furthemore, 
in [26], the regularity of solutions except for a countable set of singularities whose projection on 
the time axis is a finite set had been obtained. Very recently, in dimension three, Lin and Wang 
[28] have proved the existence of global weak solutions under the assumption that do £ §+ by 
developing some new compactness arguments, here is the upper hemisphere. 

As for the compressible case, Huang et al. m proved the local existence of unique strong 
solution of m provided that the initial data po , u o, do are sufficiently regular and satisfy a 
natural compatibility condition. And a criterion for possible breakdown of such a local strong 
solution at finite time was given in terms of blow up of the quantities |p||and ||Vd|| L 3 L oo. 
In |L8], an alternative blow-up criterion was derived in terms of the temporal integral of both 
the L°°-norm of the deformation tensor T >u and the square of the L°°-norm of Vd. In terms of 
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the global well-posedness, results in one dimensional space have been obtained in mm- In two 
dimensions, Jiang et al. [2T] established the global existence of weak solutions under the small 
initial energy. In two or three dimensions, if some component of initial direction field is small, 
Jiang et al. m established the global existence of weak solutions to the initial-boundary problem 
with large initial energy and without any smallness condition on the initial density and velocity. 
Recently, Lin et al. [29] established the existence of global weak solutions in three-dimensional 
space, provided the initial orientational director field do lies in the hemisphere §■]_. 

The low Mach number limit of the system (II . 1 D - (II .21) has also been studied recently. Hao 
and Liu [15] investigated the so-called incompressible limit (i.e., the low Mach number limit) 
for solutions in the whole space W N (N = 2,3) and a bounded domain of M. n (N = 2,3) with 
Dirichlet boundary conditions. Ding et al. [8] studied the incompressible limit with periodic 
boundary conditions in M. n (N = 2,3). Wang and Yu [35] proved the incompressible limit for 
weak solutions in a bounded domain. For more about the incompressible limit problem, one can 
refer to [221 EH E2] and the references therein. 

Let us mention that all of the above results were performed in the framework of Sobolev 
spaces. Inspired by [5] for the compressible Navier-Stokes equations, it is natural to study the 
system (PJ-CuD in critical Besov spaces. We observe that the system m is invariant by the 
transformation 

p = p(l 2 t 1 lx), u = lu(l 2 t, lx), d = d(l 2 t,lx) (1.3) 

up to a change of the pressure law P = l 2 P. A critical space is a space in which the norm is 
invariant under the scaling 


(e,f,g)(x) = (e(lx),li(lx),g(lx)). 


Very recently, in the case N = 3, Hu and Wu m studied the global strong solution to (ll.ll) - (11.2D 
in critical Besov spaces provided that the initial datum is close to an equilibrium state (1,0, d) 
with a constant vector d £ S 2 . More precisely, there exist two positive constants ijo and Tq such 


that if po — 1 £ 


, u 0 £ B 2 2 ,i(®" ), and d 0 — d £ Bj 


satisfy 


11/90 - 111-1,00 + IKII 1 T 11 d 0 d|| i < 770, 


B4 


(1.4) 


then system (ll.ll) - (ll.2l) has a unique global strong solution (p, u, d) with 


p - 1 € L 1 (R+-, B3’ 1 ) n C(R+; R|’~), 


+ • 


u€ L i (M+;B 2 MnC(M+;B 2 M , 


d — d € [ L 


\ 1 nn>+ ; rJ’ 00 ) n b|’°°) 


(1.5) 
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satifying 


IIP-1 


L°°(K 

+ IIP-1I 


-,bI~) + l|u|| i. 


i + lid — dll 

■Md L° 




( 1 . 6 ) 


3 ! + U 

L!(R+;S^ ) L° 


5 + lid — dll _ 7 „ < r 0 r? 0 . 

-;B 2 \) " L°°(R+;B^’ ) “ 


Here B* r and Bf, 1 ' 1 denote the homogeneous Besov space and hybrid Besov space, respectively. 
We are going to explain these notations in Section 2. 

The purpose of this paper includes the following two aspects: 

On one hand, we establish global strong solutions to the Cauchy problem of (ll.ll) - (ll.2l) in 
critical Besov spaces with initial data close to a stable equilibrium. From nzi, when N = 3 
and (|1.4I) holds true, the system (II. ![) - (!!.2D has a unique global strong solution (p — 1, u, d — d) 
satisfying (|1.6I) . Concerning the global well-posedness with respect to d, we carry out in the 

.3 _ 3 ]OC 

framework of critical Besov space B^ j (if N = 3) but not the hybrid Besov space B£ ’ in [17j . 

-Too -1 -3 -3 

Since B$ ’ ~ B 2 2 1 n 1 C B,^ { , the regularity requirement of the initial data in terms of the 

director field do — d is relaxed. The key point is that, different from the estimate of d in hybrid 
Besov space in m Proposition 4.1], we make use of the estimate of d in the homogeneous Besov 

. 3 

space B 21 (see Proposition 13.21 below when N = 3). In addition, the global estimates of a linear 
hyperbolic-parabolic system given by Danchin [5] (see also [U Chapter 10]) play an important 
role. 

On the other hand, we give the rigorous justification of the convergence of the incompressible 
limit for global strong solutions to the compressible equations of liquid crystals when the initial 
data are ill prepared and small in a critical space. Meanwhile, the accurate converge rates are 
obtained. Our proof follows the ideas of Danchin [6j and the key point is to use some dispersive 
inequalities for the wave equation: the so-called Strichartz estimates (see e.g., [HEHliSlj and 
the references therein). 

We would like to point out that [7] is the first paper devoted to the incompressible limit 
problem where Strichartz estimates have been used. In the spirit of [7j, Danchin [6] studied the 
zero Mach number limit in critical spaces for barotropic compressible Navier-Stokes equations. 
Fang and Zi m investigated the incompressible limit of Oldroyd-B fluids in the whole space. 

Before presenting the main statements of this paper, we introduce the following function 
space: 


K(T) = {(e, f, g) € ^(O^T^B^nC^Tj-Btn) 

/ • \ N 

x (^(O^jHl+^nCaO,^;^ 1 )) 
(^(O^;^ 2 ) nC'([0,T];5| il )^ JV | 


x 
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and 


ll( e >f’g)ll'BJ(T) — ll e ll + Ill’ll L^{B s 2 ~-f) + llgllL5?(S| :1 ) 

+ v W e \\L) r {Bt' 1 ) +^ll f lli,i,(i? 2 s + 1 ) + 6 'IIsIIli,(^+ 2 )- 

Here T > 0, s E M, u := A + 2p and v_ := millet, A + 2 p). We use the notation B s u if T = +oo 
by changing the interval [0, T] into [0, oo) in the definition above. 

Our first result of this paper reads as follows. 


Theorem 1.1. Let d E 


be an arbitrary constant unit vector, and assume that P'( 1) = 1. 


There exist two positive constants r] and F such that if po — 1 € RJ 

N_ 

Bf x satisfy 


,41-1 


uq € B 2 \ and do — d E 


IlPo - l|L#,oo + ll u o! 


+ l|d 0 - d| 


B. 


2,1 


B. 


< V, 


(1.7) 


2,1 


then the following results hold true: 

N_ 

(i) System (11.1D - (11.2j) has a global strong solution (p, u, d) with {p — l,u, d — d) in BJ 
satisfying 


ll(P — 1, u, d — d)|| n 


< r 


Po 


- II 


-4 l00 

b7 


u o| 


N 

JDr> 1 


+ l|d 0 - d| 


N 

b" 1 

JDr, -I 


( 1 . 8 ) 


(ii) Uniqueness holds in BJ if N > 3 


and in Bl n Bf, (1 < s < 2) if N = 2. 


Remark 1.1. Since P(p) is an increasing convex function of p, we assume P ; (l) = 1 for 
simplicity. The general barotropic case P'{ 1) > 0 can be verified by a slight modification of the 
argument below. 


Recall that the Mach number for the compressible flow (11.11) is defined as: 


Vnp) 

Thus, letting M approach zero, we hope that p, d keep a typical size 1, and u is of order e, 
where e € (0,1) is a small parameter. As in [30] . we scale p, u and d in the following way: 

p = p e (et,x), u = eiT(et,x), d = d e (ei, x), 


and we take the viscosity coefficients as: 


p = ep e , A = eA £ , £ = e 2 f e , 8 = e6 £ . 
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Under this scaling, system (ED-dSD becomes 


' d t p e + div(/Au e ) = 0, 

dt(p e u e ) + div(p e u e 0 u e ) — p e Au £ — (p e + A e )Vdivu e + - 

e z 

< = -£ £ div(Vd e O Vd e - ±|Vd e | 2 l), ( 1 . 9 ) 

d t d e + u e • Vd e = 0 e (Ad e + | Vd e | 2 d e ), 

k (p £ j u > d e )|t=o = ( pq , Uq, dg). 


For the simplicity of notations and presentation, we shall assume that //,A e ,£ e and 9 e are 
constants, independent of e, and still denote them as p, A, £ and 9 with an abuse of notations. 
Formally, we get by letting e —)• 0 the following incompressible model 

dtu + u • Vu + V 7r = pAu — £div(Vd 0 Vd), 
d t d + u • Vd = 0(Ad + |Vd| 2 d), 

< ( 1 . 10 ) 
divu = 0, 

(u, d)|t =0 = (u 0 , d 0 ). 


Thus, roughly speaking, it is also reasonable to expect from the mathematical point of view that 
the global strong solutions to ( 11.911 converge in suitable functional spaces to the global strong 
solutions of (11.101) as e —>• 0, and the hydrostatic pressure ir in the first equation of (11.1011 is 
the limit of ^ — 4-|Vd e | 2 . Our second goal is devoted to the rigorous justification of the 
convergence of the above incompressible limit in the whole space. We remark that the existence 
of global strong solutions to the incompressible flow of liquid crystals ( 11 . 101 ) in critical Besov 
space was established in Xu et al. [ 35 ] - 

As in [6], we want to consider so-called ill prepared data of the form p e Q = 1 + eb g, Ug and dg, 
where (i>g, Ug, dg) are bounded in a sense that will be specified later on. Setting p 6 = 1 + eb e , it 
is easy to check that (6 e ,iT,d e ) satisfies 


diviT 

8 t b e + = —div(6 e iT), 

_ - f f pAu e + (p + A)Vdivu e P'(l + eb € ) V6 e 

1 + eb e 1 + eb e e 

= div ( Vd e © Vd e - -iVdd 2 !) , 

l + eb e V 2 1 /’ 

d t d € + u e ■ Vd e = 0(Ad e + |Vd e | 2 d e ), 

. (6 € , W, d £ )|t =0 = (6g,Ug,dg). 


( 1 . 11 ) 


Our second result of the paper can be stated as follows. 
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Theorem 1.2. Denote Q := VA _1 div and V := I — Q. Let d E be an arbitrary constant 

~K,00 . *L-1 

unit vector. There exist two positive constants c and M such that if £ By ' , Ug E 
~ . K 

and dg — d E Bf x satisfy (for all 0 < e < eo) 

( 1 . 12 ) 


■— H^oII -4-i + ei/ ll&oll . 4 + ll u oll. 4-1 + Ho d11 4 < c, 
°2,1 °2,1 -°2,1 ^2,1 


£/ien the following results hold: 

1. Existence: 

N_ 

• System ([1,111) has a solution in 03 e i such that for all 0 < e < eo, 

||(6 e ,u e ,d e - d)|| n < MCq p . 

m2 

^eu 

• System (11.101) has a unique solution such that 


u 


n ,, “f* 11 d — d|| n T 11 d — d || n , 2 

'L°°{b7i *) " l-HflaTi ) 


.#-1 + ll u ll 


< M (l|u 0 || 4_! + ||do — d|| N 
^2.1 ^2,1 


2. Convergence: 

• If N > 4: For all p E [pat, oo] with pn := 2( ^_ 3 1) , we /iare 


. ,A_1 + II Qu C II _ . JV 1 < MCg^eS , 

2 ) 2 ) 


and 


||Fu e — u11 jv_3 + ||Fu e — u|| jv i i + ||d e — d|| 

L°°(B? 2 ) LHB p L°° (Bf, 2 ) 


N__ 1 


5 ) 


+ ||d e -d|| . jv + s < Ml ||Fuq — u 0 || _ jv_3 + ||dg — d 0 || _ jvi + Cq u 


B. 


jr~? 

p, i 


B. 


p, i 


If N = 3: For all p E [2, oo), we have 


2 p 21 +||Qu C || 41 < MCq^C 2 P, 

L*=*{BZ i 2 ) i 2 (^i 2 ) 


and 


||Fu e — u|| 4 _ 3 + ||Fu e — u|| 4,i + ||d e — d|| 4 _ i 


+ ||d c -d|| 4 + 3 < M( ||Fuq - u 0 1| 4_3 + Hq - d 0 || 4_1 + Cq u e 2 

lhb;,*) v bl^ 


i_i 

p 


B p 72 

n p,l 


If N = 2; For all p G [2, 6], we have 


4 2 _ . 4 —3 + ||Qu e ||_ ^-l < MCq V €* i ; 


*) " "LHb% 4 ) 
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and 


||’Pu e — u1 1 5 _5 + \\Pxi e — ull _ __ 

L°°(B^i *) Li(B^ *) 


5_, 3 + ||d e — d| 

\ 

P, 1 


*) 


+ " de_d| WA +| ) - M( " Pu °“ u oi 


A_5 + || dp — do 


P,1 


I s i+Cfe^ 
° 

p, 1 


Remark 1.2. When d = d, the compressible flow of nematic liquid crystals is reduced to the 
well-known compressible Navier-Stokes system. Our results coincide with the ones in [6| con¬ 
cerning the zero Mach number limit in critical spaces for compressible Navier-Stokes equations. 


Our paper is organized as follows. In the next section, we recall some basic facts about 
Littlewood-Paley decomposition and the homogeneous Besov spaces. In Section 3, we investigate 
the existence of global solutions for system (HMD-(H3D- Section 4 is devoted to the proof of 
uniqueness. In Section 5, we will prove the convergence of the incompressible limit in the whole 
space . 

We end this section by introducing the notations used throughout this paper. C stands for a 
harmless constant which never depends on e, and we sometimes use the notation A < B as an 
equivalent to A < CB. The notation A~ B means that A < B and B < A. 


2. Homogeneous and hybrid Besov spaces 

We first recall the definition and some basic properties of homogeneous Besov spaces. They 
could be defined through the use of a dyadic partition of unity in Fourier variables called ho¬ 
mogeneous Littlewood-Paley decomposition. To this end, choose a radial function ip € S(M . N ) 
supported in C = {£ £ M. N , | < |£| < |} such that 

E <p(2~ j 0 = 1 if Z* 0. 

je z 

The homogeneous frequency localization operator A j and Sj are defined by 

A ju = D)u, SjU = AfcU for j € Z. 

k<j-l 

With our choice of one can easily verify that 

A p A q u = 0 if \p — q\ > 2 and A p (S q _iuA q u) = 0 if \p — q\ > 5. (2-1) 

Let us denote the space /(l^) by the quotient space of S'(K N )/V with the polynomials 
space V. The formal equality 

U = Y1 A kU 
fcez 

holds true for u € and is called the homogeneous Littlewood-Paley decomposition. 

We will repeatedly use the following Bernstein’s inequality: 


GLOBAL EXISTENCE AND INCOMPRESSIBLE LIMIT FOR FLOW OF LIQUID CRYSTALS 


9 


Lemma 2.1. (see 0) Let C be an annulus and B a ball, 1 < p < q < +oo. Assume that 
f e LP(R n ), then for any nonnegative integer k, there exists constant C independent of f, k 
such that 

supp/ C A B => \\D k f\\ Lq(MN) := sup \\d a f\\ Lq( ^ N) < C k+1 \ k+Ni p~* ] \\f\\ LP{R N ) , 

\a\=k 

supp/ C AC ^ C-*- 1 A*||/|| l , (r * ) < \\D k f\\ LP{RN) < C k+1 X k \\f\\ LP{RN) . 

Next, let us recall the definitions of the Besov spaces. 

Definition 1. Let s € M, 1 < p, r < +oo. The homogeneous Besov space B p , r is defined by 

B s p ,r = {/ € y\R N ) : ||/||^ < +oo}, 

where 

ll/ll B- :=l|2 fa ||A t /||„|| P . 

P, T 

We next introduce the Besov-Chemin-Lerner space L^(Bf r ), which is initiated in [2]. 
Definition 2. Let s € R, 1 < p, p, r < +oo, 0 < T < +oo. The space Lj,(B^ r ) is defined by 
L P T (B s p ,r ) = {/ € ( 0 , +oo) x y'(R N ) : ||/|| z , ( ^) < +oo}, 

where 

ll/ll L p T (B^ r ) := 1l Z ^fc/(i)||LP(0,T;LP)||r- 

A direct application of Minkowski’s inequality implies that 

L p T (B* Ptr ) ^ 4( B px). ^ r>p, 

L p T {B s p , r ) ^ L p T (B s p J, if p>r. 

We also need the following hybrid Besov space introduced by Danchin in [5]: 

Definition 3. For u > 0, r£ [1, +oo], s € R, we define 

3 6Z 

By the definition, it is easy to verify that 



which means that Bf’°° ~ i 1 Cl 

Let us now state some classical properties for the Besov spaces. 
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Proposition 2.1. The following properties hold true: 

1) Derivation: There exists a universal constant C such that 


C 


-i 


B 3 


< IIV/IU-l < C\\f\\g. 


2) Sobolev embedding: If 1 < pi < P 2 < oc and 1 < r\ < r<i < oo ; then 


7*1 


N_ , N 
Pi 
P2,r 2 


TD S Pi ' P2 


3) Real interpolation: ||/|| ■ e sl+ (i_e) S2 < 


R S 1 

J-'p^r 


l-e 

R 3 2 ' 
^p, r 


4) Algebraic properties: for s > 0, B* ) X n L°° is an algebra. 

5) Scaling properties: 

(a) for all A > 0 and f S Bp^, we have 


ll/(A-)llj3.*A-7 

p, i 


B s . > 

p,i 




(b) /or / = f(t, x ) m L r (0, T; !//,), we /iare 

l|/(A“-,A 6 -)||_ L r(B» jl ) « A - ' IU 

Next we recall a few nonlinear estimates in Besov spaces which may be obtained by means 
of paradifferential calculus. Firstly introduced by J.-M. Bony in [2], the paraproduct between / 
and g is defined by 

Tf9 = ^ ^Sq-ifAgg, 


and the remainder is given by 


with 


R{f,g) = '52& q f\g 


^qd (^g-1 + A q + A q+ \)g. 
We have the following so-called Bony’s decomposition: 

fg = T g f + t f g + R(f,g). 


( 2 . 2 ) 


T' 


f9 


The paraproduct T and the remainder R operators satisfy the following continuous properties. 
Proposition 2.2. Suppose that s € R, a > 0, and 1 < p,pi,p 2 ,r,r\,r 2 < oo. Then we have 


1) The paraproduct T is a bilinear, continuous operator from L°° x Bf r to Bf r , and from 


B, 


00,7*1 


x B s p r2 to B s v / with A = min{l, ± + £}. 


P,T 


7*1 7*2 - 


2) The remainder R is bilinear continuous from Bf l i ri x B££ r2 to Bf l + S2 with si + S 2 > 0, 
= — + — < 1, and - = — + — < 1. 

From (12.21) and Proposition 12.21 we have the following more accurate product estimate: 
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Corollary 2.1. If u E B p * , and v E B p 2 2 , with 1 < p\ < P 2 < oo, si < S2 < ^ and 

■ S\+S 2 — — 

si + «2 > 0, then uv E B p , t 1 P1 and there exists a constant C, depending only on N, s i, S 2 ,pi 
and p 2 , such that 


\\uv 


I +a2 _iL < C\\u\ 

b p 2 a 


R S 1 

B piA 




R s 2 ■ 
B p 2 A 


We finally need the following two composition lemmas (see mm)- 


Lemma 2.2. Let s > 0, p E [1, oo] and u E B pl n L°°. Let F E l / F 1 ^ i ! +2 ’ 00 (M Ar ) such that 
F( 0) = 0. Then F(u) E B s p l and there exists a constant C = C(s,p, N, F, ||ii||x,oo) such that 


IIWIIb. <c\\u 


p, i 


p.i 


• — • riV] , o 

Lemma 2.3. If u and v belong to Bf l; (v — u) E -Bf,i /of s E (—^ , t! anc ^ ^ € ^/oc ' 0^) 
satisfies G'{ 0) = 0, then G(v) — G(u ) belongs to .Bf l an d there exists a function of two variables 
C depending only on s, N and G, and such that 


II G(v)-G(u) 


I j ^ (IMI-L 00 i 


\V\\L°° 


N 

r"2" 

■ d 2,1 


+ M 


JV 

r"2" 

i 3 2,1 


n — u 


B: 


2.1 


3. Global existence for initial data near equilibrium 

In this section, we will prove the part (i) of Theorem ll.il Throughout this paper, we consider 
only viscous fluids, those for which p > 0 and v > 0. The following proposition plays an 
important role in obtaining the estimates of (u, b). 


Proposition 3.1. (see ID) Let Vv E Lj,(L°° (M N )), s E M, and (a, u) 6e a solution of the 
following system 


dta + div(T v a) + divu = F. 

dtu + T v ■ Vu — ^Au — (p + A)Vdivu + Va = G, 


(3.1) 


on [0, T). Then there exists a constant C depending only on N and s, such that the following 
estimate holds on [0, T): 


bF 


> + 


l u (*)ll Ki 1+ J 0 ( 


is\\a\\ $3,1 + v\\ u \\ b s+1 


L )dt' < C' e C||Vv|l G 1 ^°°(K JV )) 


x ( ll^ollo?.. 00 + IIuoIIb-t 1 


2,1 


f 


‘ e - C " Vv l t; , (t » ( R «„ (||F|| | |G|| )dt , 


(3.2) 


Next, we establish the estimate of the director field d in critical Besov sapce. 


Proposition 3.2. Let s E (— y, 1 + ff), and d be a solution of the following equation 


dtd + u • Vd — 6>Ad = M 


(3.3) 
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* —+1 

on [0, T), where u € Ll r (B 2 I ] ) and M. G Ly(S| i). Then there exists a constant C depending 
on N and s, such that the following estimate holds on [0, T): 


i) + ^ll c ^llz/i(s|+ 2 ) — C* ex P 


U|| N,, 


( 3 - 4 ) 


Proof. Applying A q to (13.31) and taking the L 2 inner product of the resulting equation with 
A q d, integrating by part, we deduce that 

~||A 5 d||^ + 9||vA,d||^ 

1 


< ^-||divu||i<x>||Aqdll^ + ||[u, A,] ■ Vd|| L 2 + \\A q M\\ L 2 )\\A q d\\ L 2. 
Hence, according to Bernstein’s inequality, we get, for some universal constant k, 
\\\^\\ l ^°( l 2 ) + ^ K 2 29 ||A g d|| L i( L 2) 

< ||A 9 d 0 || L 2 + J ^||divu||L<x>||A 9 d|| L 2 + ||[u, A q ] • Vd|| L 2 + ||A q M\\ L 2^dt'. 
Now, multiplying both sides by 2 qs and summing over q, we end up with 

HIIlHbIu) + Ke \\ d h}(B s 2 + 2 ) 

< lldolls^ + (^l|divu||i<x;||d||^s^ ^ 2 « s ||[u, A ? ] ■ Vd|| L 2 + WMW^ 


dt!. 


(3.5) 


According to Lemma 2.100 in [I], we get, for— y < s < ^ + 1 j 

^2^||[u,A 9 ].Vd|| L 2<C||u|| Y + 1 ||d|| B 

-^2,1 




(3.6) 


Substituting (13.61) into (13.51) and using the embedding Bfi L°°, we get 

H d llit°(4 s ,i) + Ke W d hl(B s + 2 ) 


- + C j 0 


a+i ||d ||gs dt, i 


b" 1 

n 2,l 


2,1 


taking advantage of Gronwall’s inequality, we obtain (13.41) immediately. This competes the proof 
of Proposition 13.21 □ 


Proof of the part (i) of Theorem 11.11 We are going to prove that if the initial data 


IIpo - !|| - 4,00 + ll u oll -4-1 + ll d o - d|| N < 77 , 


B. 


Bo 


(3.7) 


for some sufficiently small rj, there exists a positive constant T such that 

||(p- l,u,d-d)|| n < Trj. 

537 


(3.8) 
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This uniform estimates will enable us to extend the local solution (p — 1, u, d) obtained by using 
a Friedrichs method as in [1] to be a global one. To this end, we use a contradiction argument. 
Define 

T 0 = sup |t £ [0, oo) : ||(p - l,u,d- d)||^jv < Ft] j (3.9) 

with T to be determined later. 

We note that Vd = V(d — d) because d £ is a constant vector. Letting b = p — 1, we 

rewrite the system m as 

dtb + divu + div(T u 6) = —div(T t / u), 

< d t u + T u -Wu-Au + Vb = K(b)Vb-I(b)Au-f! )iU u i + H, (3.10) 

k d t ( d - d) + u • V(d - d) - #A(d - d) = J, 

where 

H := - T ^_div(v(d - d) 0 V(d - d) - i|V(d - d)| 2 l), 


J := 0|V(d - d)| 2 (d - d) + 0|V(d - d)| 2 d, 

and A := pA + (p, + A)Vdiv, 1(b) := K(b) := 1 - P . 

Suppose that To < oo. We apply the linear estimates in Propositions 13.11 and 13.21 to the 
solutions of reformulated system (13.1011 such that for all t £ [0,To], the following estimates hold: 

\\ b (t)\\ -4,00 + ll u WII .f-1 + [ (v\\ b \\ -f.1 +lI|u|| N +1 )dt' 

B .7 B7 n v BJ b7 / 


< 


Ce "4 0 (iOO(RjV)) {||5 0 || ^ jOQ + ||uo|| n_ x + ||div(T b 'u 


l'(b7’ ) 


- —1 1 II- 1 jv T ||Tq. 1] U*|| N _ 

LUBZ ) ; ) 1 Lh (B7 *) " *) 


+ \\Kmb\\ » + \mAuu 

^T 0 (^2,l ) 


1 T ||iL|| . 4v._ t 1, 


(3.11) 


and 


l|d(t)-d|| n +0 

T2 o i , 


fT 0 


Id — dll N, 0 dt' 

b t+ 

-°2,1 


< Cexp ( C||u|| m 

2,1 ^ 


do — dll at + 

r“2" 

■ d 2,1 


(3.12) 


. JV 

^T n (^20) 


J 0 ' A* J o' 

In what follows, we derive estimates for the nonlinear terms one by one. Similar to the case 
of isentropic Navier-Stokes equations pTj, by Proposition 12.21 Corollary 12.11 and Lemma 12.21 we 
have the following inequalities: 


l|div(rWll u^L", 5 


N U N ,, , 

l-(r7’ ) ) 


(3.13) 
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||A'(6)V6|| 


L —i ^ 




||/(6)^ U || N, 


< 


\\T' di y 


-i, Z ll u ll 


N „ U 


N, U 

T h 11 1 



(3.14) 


(3.15) 

. N +1 • 
l t 0 +24 ) 

(3.16) 


Compared with the isentropic Navier-Stokes equations, the new terms can be estimated as 
follows: 


\H\ 


N , = 


^t 0 ) 


1 + 6 


div 


(v(d - d) © V(d - d) - V(d - d)| 2 I 


^T n (^2“l ) 




div^V(d - d) © V(d - d) - i|V(d - d)| 2 I 


. N 

^T 0 (+24 ) 


+e 

^ (i + 
^ (i + 

< fi + 


1(b) div(V(d - d) 0 V(d - d) - ^|V(d - d)| 2 l) 


2U ) 


(3.17) 


N 

££+ 24 ) 


V(d - d) © V(d - d) - -IV(d - d)| 2 I 


. N 

BT n +24) 


- N 


)||V(d-d)|| 


^t 0 +24 ) 


Id — d|| n ||d — d|| at , 

+-(b7’ )/ A-(43I) ) 


, = d|||V(d-d)| 2 (d-d)|| , + 0|||V(d-d)| 2 d|| * 

■^ Tq '•■^' 2,1 / ^ Tq v -^ 2,1 / ^ Tq *+* 2,1 J 

For the estimate of J\, it follows that 

J 1 =0|||V(d-d)| 2 (d-d)|| N 

L Tq 


Ji + J 2 . 


(3.18) 


< 

r^j 


|||V(d — d)| 2 1| . iv || d — d 11 

l t 0 (+ 24 ) B%? o (+ 3 ) 


N 

ill 


< ||d — d11 N ||d — d11 n 

L 2 r 0 (BZ ) L ¥ 0 ( B &) 


(3.19) 


<||d —d|| . 4+2. I|d — d|| n 

L T n ( B 2,l ) L T 0 (B^!) 


'O' 


For J- 2 , we have by the definition of Besov’s spaces 

|2 -ll 

+0 


J-2 = <9||IV(d — d)| 2 d|| n < ||d-d|| 2 

L Tn (B 2tl ) L T n (B7i ) 


AT , 1 

t+g 


< lid-d|| 


B¥ 0 (b£) 


n lid — d|| 

L T 0 (B?r) 


4+2 x - 


(3.20) 
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Plugging inequalities (I3.13I) - (I3.20I) in (13.111) and (|3.12l) . we thus get 


and 


l™(b7' ) 


+ u 


L-(B; 


< 


To'-' ' ~ 1 o 

H^ u llz,l, (l°°(r n )) 


. 4-1 + [ ( V M - 4.1 + tIMI - 4 +1 W 
) J 0 v b7’ R Z J 


b„ 


. 4 ,°° + ll u o|| . 4-1 + 


So 


JV „ U J/ M 

) £k(B7i ) 


+ 


Ji + U 


iv „ 0 _ i\ , ...... . 


|U|| jv 

r 0 


+ 1 . 


+ 1 + 


L-CRj’ 00 ) 


Id — d| 


l t 0 ( r Zi) 


v nd — dll 


iV , o 
?- 2 -+ 2 ' 


}• 


Id — d| 


N 


/•To 

,+e 

1 j 0 


|d-d| 


4.4-2 

rW+ 2 

"2.1 




^ exp 


u 


^ 0 (4Y) 


){' 


do — dll iv + lid — dll iv 10 lid — dl 

LUbZ + T 




+ II d — dl 


n ,0 d — d jv 

^(B 2 Ti ) l tZ b Zi) 


}• 


The above two inequalities combined with the definition of To yield that 
II (6, u, d-d) II JV < Cie Cir?? (T7 + ry+ rV). 

*7 (To) 


We choose T = AC\ and rj > 0 satisfying 

„C\Yti ^ o t^2 


< 2, Y z r) < TV < i 


Consequently, it follows from (13.231) and the above choices of T and 77 that 


(3.21) 


(3.22) 


(3.23) 


(3.24) 


||( 6 , u, d — d) || n < Yr), 

(To) 

which is a contradiction with the definition of Tq. Hence, we can deduce that To = 00 . Global 
existence is thus proved. □ 


4. Uniqueness 

In this section, we will prove the uniqueness of the solution to system (HD, i-e-, the part (ii) 
of Theorem mo 

Proof of the part (ii) of Theorem 11.11 

JV 

Case N > 3: Assume that ( 6 j, u*, dj)j=i 2 hi solve (11.11) with the same initial data. 
Denote 


5b = &2 ~ bi, 5u = U 2 — ui, dd = do — di. 
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Then (56, 5u, 5d) E *B „ 2 (T) for all T > 0. On the other hand, since (bi, Uj, di)i=i i2 are solutions 
to system (11.111 with the same initial data, (5b, 5u, 5d) solves 

dtSb + div(T U2 • 5b) + div5u = 5F, 

< d t 5u + T U2 -V5u-A5u + V5b = 5G, (4.1) 

<9i5d + U 2 • V5d — #A5d = 5M., 

with 

5F = — <5u • V 6 i — 6idiv5u — dh^T^t^), 

5G = -T' d . 5u u 2 — 5u • Vui — ( I(b 2 ) — I(bi))Au 2 — I(bi)A5u 


+ {(Ifa ) - /(&i))div(v(d 2 - d) © V(d 2 - d) - i|V(d 2 - d)| 2 l) 

+ —2— divfv(di - d) © V5d + V(d 2 - d) © V5d 

1 + &i v 

- il(V(di - d) : V5d + V(d 2 - d) : V5d)) + V(/C(6 2 ) - /C(6i)), 


5.M = -5u- Vd x + fl(jV(di -d)| 2 5d + (V(di - d) : V5d + V(d 2 - d) : V5d)d 2 ), 

where we used the notations A : B = l and /C(z) = / 0 ~ K(y)dy. Applying Proposi¬ 

tions EH and 13.21 to the system D, we have 


||(56, 5u, <5d)|| jv 

(T) 


< c exp ( C||u 2 


' l t(4 T i +1 ) 




+ ||5G|| ^_ 2 + ||5A4|| ) 

L \,(BJ 2 ) V 


p 


(4.2) 


. I . AT , . TV 

As in [5], we could get that 6 * E C(M + ; B 21 ) (i = 1, 2). Indeed, 6 o E B 2 2 i , bi E L°°(R + ; I ? 2 2 1 ) 

. jv, ■ JY_i • HL 

and df b, = —u, • V 6 , — divuj — 6 jdivu, E L 2 (R + ; S 2 2 1 ). And because bi E C(B 21 ) n L™ C B 21 , 

we have bi E C ([0, +oo) x M. N ). On the other hand, if (11.71) is satisfied for some rj suitably small, 
we have 

\bi(t,x)\ < - for all t > 0 and x E M N . (4-3) 

In the following, we estimate the terms 5F, 5G and 5M. respectively. Firstly, according to 
Proposition 12.21 and Corollary 12.11 one gets 


INII, 


N 

- tt— l,oo nsj 


< ||5u| 


LUBh 




+ INI 


L\,(Bj ) 1 ^T‘ K ‘ J 2,l) 

Similarly, the terms of 5G could be estimated as follows: 


L™(B, 


4 -1 ’ oo ' 1 


l u 2| 


lUb, 


N ii 

* a Tl + ) 


(4.4) 


\^diSu u 2\ 


N _ 

g t {bz ;') 


-2. iS II^HI 


N _o U 2 _ 

L 2 t {B^ ) L*{B T) 


< 


IIH 


A-u 11 u 2 


; ) l 2 t (b£) 


(4.5) 
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|| 5 u • Vui|| jv_ 2 < ||HI 4-1 IIVU! || < HHI N_ ||ui|| N ( 4 . 6 ) 

l t (b ) L t {B 21 ) L T (B 21 ) L t (B 21 ) L t (B 2 .i) 


(I(b 2 ) - I(bi))Au 2 N _ 2 < 11/(62) - /(6i)|| 4-1 Mu 2 || .f-i 


~ (^1 + 
~ f 1 + 


N + ||6 2 || JV J || V - U 2 II at _! ||( 56 || JV_J 

l t( b £ l) L5?(B^i)/ ) 


. 4 + H62II .4 ) ||u 2 1| . jv +1 ||< 56 || 


- _l_i 11 11 AT _-I < 

) i§p(s a yi ) 


( 4 . 7 ) 


l|/(6l) ^ u| Wr 2 )" 


.4 J|V 2 5 u|| <||&i|| f J|6u 

L-(B 2 ^) ^(BaTi ) TocmT' 


I ,1 JV 

L™(B^) L^Bj i) 


( 4 . 8 ) 


?(/(6 2 ) - /(6 1 ))div(v(d 2 - d) © V(d 2 - d) - ^|V(d 2 - d)| 2 l) 


^(B 2 Ti ) 


< 

r^i 


\\I(b 2 ) ~I(bi)\\ ,n_ x |||V(d 2 -d)| 

L t( B ^1 ) 


N 

l t( b ^i) 


( 4 . 9 ) 


y (1 + 


N +H 62 H 

L T ( B ^l) 


N 

l t( b ^: l) 


do — d| 


INI 


,1 11 — 11 . N , 

24 ) B t( B ^1 ) 


1 + 61 


div 


(v(di - d) © V6d + V(d 2 - d) © V6d 


2 ) 


< 

nsj 


< 

nsj 


! + ||6iIL„,,4. )||V 5 d|| .4^ (||V(d x -d)|| 


Lt&Ziv l K b Z i ) v " ' * '" B U B Z I) 




v 2,1 / 

Similar to ( 14 . 101 ) . there holds that 


1 + HM n ) ||< 5 d|| N ( ||d! — d|| 

11 L h ¥ (B^)J" "l^b^v" 4(®2 T i T1 ) " “ "yyn 


jv, 1 


^4 +||V(d 2 -d)||^^N.J ( 4 . 10 ) 

l) 

d 2 - d|| 


N, 1 


—divf-I(V(di - d) : V 5 d + V(d 2 - d) : V 5 d 
1 + b\ V 2 


ii.(B a Ti ) 


(1 + 


1 V 1 1|oai| JV Mid!—dll JV,, -I- 110-2 — an JV,, 

l t( b Z)' b t(b 7 i)^ ^(B 2 V ) B?-(B a y i + ) 


+ ||d 2 - d|| 


(4.11) 


Moreover, applying Corollary 12.11 and Lemma 12.31 one easily gets 


IIV(/C(6 2 ) - /C(6!))|| 4 < ||/C(6 2 ) - K(h )|| 4 

L t \B 2,1 ) ^tI° 2,1 ) 
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Note that in the above inequalities, we have used that N > 2. Collecting the above estimates 
(14.5f) - (14. 12[) . it follows that 


\\SG\\ 


N _ o 

LH B ^1 ) 


< I ||Ul|| JV + 11U211 JV I lldull N_ 1 + 


JV lldull JV 

l) l t( b Z) 


+ l 1 + 

+ ( 1 + 

+ (i + 

+ T 


N 


+ INL • t s ) II U 2 


- - JV +1 j|56|| JV _J 

l t ( s 2fi) L ¥(B^i)J ) L T (®2U ) 


JV 


.4 + IIM 


n I ll u 2 

l t ( B 2u) 


d 2 -d|| 2 jv +1 ll^ll . JV 

) L t(^1 ) 


JV 

l t ( b Zi) 


di — d| 


L UB? : L +1 ) 


+ ||d 2 - d| 




v 2,1 


+ IIM, 


INI, 


JV 


' L T (^2u) L T C^l) 7 " " L T ( B 2 2 l d 

For the estimate of &M, note that 




(4.13) 


||5u - Vdi|| = ||<5u-V(d x -d)|| .n^ 

LU B Z ) L t( B Z ; ) 


< 


IIHI ||V(di -d)|| n < IIH 

M4,1 ) b t (b 21 ) 


jv _, ai 

L 2 A B Z ) 


di — d| 


+1. ’ 


'lUbZ^) 


(4.14) 


IIW* - d) l a “ W,-, s l|M W,-> ll|v(dl - d)|2|l n<B?.» 


< 


||5d|| * , ||d!-d| 

L ¥( B Zl ) 


. 2V+1 > 

lU b Z ) 


(4.15) 


||0(V(di-d): V<5d + V(d 2 -d): V<Sd)d 2 || x .n_, 

< ||(V(d x - d) : Vtfd + V(d 2 - d) : V5d)(d 2 - d)|| 


N 

la b Z ) 


+ II (V(di - d) : V<5d + V(d 2 - d) : Wd)d|| . n 


lAbZ ) 


< lido — d| 


di — d| 


u 2 — an 1 v 1 11 ui u 11 jv 1 , + || d 2 — d || jv 1 , .. 

l t( b Zi)' LA B Z ) L A B Z )■’ l t( b Z ) 


(4.16) 


+ l|di-d| 


jv,, + ||d 2 — dh jv,, llioaii n 
L^(BZ + ) "r2jnT+V" "fSjqTi 


L 2 T ( B Z V L U B Z) 


< (1 + IId 2 — d| 


JV 

l t( b Zi ) 


di — d 11 jv,, + ||d 2 — d 11 jv,, 1 iiouii jv 

la b Zi ) la b Zi v L2 A B Zi) 
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Therefore, 

\\6M\\ 


N 

L t( B Zi ) 


< 


||< 5 u|| i v 1 lldi — dll jv +1 + ||< 5 d|| n j Ildi — d| 

^(4T ) ; +1 ) V(4Y) 


iV 1-1 


(4.17) 


+ 1 + || d 2 — d|| .iv 


di — d|| jv.j + ||d 2 — d|| jv,, )||M|| at . 

i|(s 2 V ) ^(S 2 V V b Hb 7 i) 


Consequently, it follows that 


||(<SMu,ad)|| * , <Z(T)\\(5b,5u,5d)\\ , 

53j (T) (T) 


with 


Z(T) = Cexp(C||u 2 || n +1 ){„ 

) J 1 l%{b7' ) 


+|U! W i , +|Ul1 


N 

l 2 A b Z) 


+ ||u 2 || N +(l+||£>l|| N + II ^2 II N 

L t( B 2l ) ' L t( B ^i) L t( B ^1 


X Hl u 2|l r ^ +Hdi-d 

v^2,l 


iv ,-+- nui 

) "LHB^r) 


n _ + i + 11 d 2 — d| 


L^ B Z +1 ) 


+ T 


n + ||6 2 || n ) + ( 1 + ||d 2 — d|| i v + ||&i| i v 

L^(B7\) ' l t( B ?i) l t( B 7 i) 


X lid, -d| 


+i. 


4( fi 2 T l ) 


+ ||d 2 - d|| 


+i, 


m B Zi ) 


It is now clear that 


limsup Z(T) < C < ll^ 1 || ron , -^,o 


T—>-0+ 


£?W’ ) 


Thus, if r) > 0 is sufficiently small, we have 


||(«5Mu,M)|| n, =0 

Bj (T) 

for certain T > 0 small enough. So (&i,ui,di) = (6 2 ,u 2 ,d 2 ) on [0,T], 

Then we could argue as in [5] for the isentropic compressible Navier-Stokes equations. Let 
T m be the largest time such that the two solutions coincide on [0 ,T m \. Taking T m as the initial 
time, we define 

( bi(t ), u i(t), d i(t)) := ( bi(t + T m ),Ui(t + T m ), d*(i + T m )). 

Repeating the above arguments and using the fact that ||6j(0)|| 

L°° < \ (see m) , we can prove 

that 


ui(i), djV)) = (6 2 (t),u 2 (t),d 2 (t)) 
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on a suitably small interval [0, e] (e > 0). This contradicts the assumption that T rn is the largest 
time such that the two solutions coincide. As a result, T m = oo, which means that the uniqueness 
result holds in M + . 

Case N = 2: Using again Propositions 13. ll and 13.21 it follows that 


\\(5b,5u, <Sd)||,g.-i (r) 


< 


Cexp (C'||u 2 || 1 u ( £ 2 i) ) (ll^llii + \\5G\\ l ^ s 2 - 2 ) + II^IIli,^- 1 ))- 


(4.18) 


Let us estimate 5F,5G and 6JA, respectively. 

Il^ll L^Sr 1 ’ 00 ) ~ T ^W 6u K^,*i Jl b illL3?(fl;-») + W^Wl^b’JMl* 

+ ll^llL-(sr 1 ’ 00 )ll U2 H^(4 2 i) 

< r^|| < 5u||^^. i) ||<5u||^ (A3 _ 2) ||6 1 || jDS?(j g-.oo ) + IIHL^) 

+ ll^llL-(Br 1 ’ 00 )ll U2 ll^(S|,i) 

< T 2 (j|<5u|| L i^ j + II^ u IIl5?(B“- 2 )) HM 

+ 11 11 ;1 ) H^l + ll^ll i; ~(BS- 1 ^)ll u 2|| z ^(B2i), 




where we have used the following inequality: 




Moreover, 


H^H^^-2 


< 

r^j 


(|I Ui IIl2 ( b ^) + H^IIlkb^)) II^IIlks- 1 ) + II 6i IIl-(b 2 i 1 )II <5u IIli,(A| ;1 ) 

+ ll 62 lli??(S 2 1 , 1 ))ll U2 11^(51,!) 

m\i 


H 1+ 

+ ( 1+ 

+ (l + 


L?(B 2 (i) + ll & 2 llL«(flJ il )jll d 2 

L ?(Bi i)) (ll d l “ d HL^i) + H d2 ~ d Hi|(Bl,i)) II^Hl 2 (B| 4 ) 
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ll^^ll 
£ P u 


^(B5,1 1 )H dl_d lli?’(B a a 1 i) + ll <5d ll 


L-(R 2 ^ 1 )H d l d H L^BU 


+ { 1 + l|d2 “ d Hi5P(B a 1 , 1 )J Ul dl ~ d|,L K4 2 ,i) + l|d2 “ d HL|(S 2 2 1 )J H <5d H^(SS, 1 )- 

Consequently, it follows that 

||(< 56 , < 5 u, ( 5 d)|| 5B .-i ( - r) < £(T)||( 5 &, 5 u, 5 d)||, 8 .-i (r) 


with 


Z(T) = Cexp ( C||u 2 || i l i( _B2 i) ^ | l|& 1 ||/oof’Rl,°o'| + l|U2||^l f R2M + ll u il 


+ II u 2|Il2 (41 :1 ) + (! + 


x ( 11 u 2 1| r,i 1 1?2 n + ||di -d|| 2 2 ^ 2 ,+ ||d 2 -d" 2 


L t ( B 2,l) 

L 2 T (Bl i) 

+ T (llM 

L 2 T (Bl i ) 


l5?(b;- oo )Tii^iili ( b 3 ii) 

OOfRl 1 + 11 M £00 (£1 A 


L2(B 2 i) 


i|(s 2 2 i) 

+ ( 1 + ||d 2 — dll roo /oi x + 


2 u Hi|?(Bj il ) I - 1,01 ULfiBlJ 


X ( ||di - d|| r a r ^a^ + ||d 2 - d|| £2 (b| l} 


+ T 


l— — 


L5?(R, 


; , 00 )}. 


It is clear that 


limsupZ(T) < CH&il^oo/^i.oow 
T—>0+ tK " 

Thus, if rj > 0 is sufficiently small, we have 

11(56, <5u,(5d)||« 8 .-i (T) = 0 

for certain T > 0 small enough. So (6i,ui,di) = (6 2 ,u 2 ,d 2 ) Q n [0, T], We can now achieve the 
proof as in the case N > 3. 


□ 


5. Incompressible limit 


This part is devoted to the proof of Theorem 11.21 Firstly we list the global well-posedness of 
incompressible model (11.101) . 

Proposition 5.1. (see }36| ) Let N > 2. Suppose that the initial data (uo,do) belong to 
B 21 (M. n ) x B 21 (S^ -1 ) with divuo = 0. Then there exist two positive constants c and C, 

depending on N,fj,,£ and 9, such that if 

||u 0 || . N ! + ||d 0 — d|| . < c, 

^2,1 ^2,1 

then there exists a unique global solution (u, d) to (|1.10|) in the class 

N_ _1 JV , i N iV i o 

^([0,00);^ )nL\[0,ooy,B 2 y ) X L°°([0, oo); n ^([0, oo); ). 






22 


Q. BIE, H. CUI, Q. WANG, AND Z.-A. YAO 


Moreover, the following estimate holds true 


u 


. + || u || 

3 ^ b r i 


n ,, “f" lid — dll n “l - lid — dll iv_i_Q 

L t( B 2 2 1 A L t( B ^1 ) L T O^fl) L t( B ?, 1 ) 


< c (||u 0 || n_ x + ||d 0 - d| 


B Z 


B Z 


(5.1) 


The study of incompressible limit problem of (fT9l) relies on Strichartz estimates for the 
following system of acoustics: 

Proposition 5.2. (see IB) Let ( b , v) be a solution of the following system of acoustics: 

dtb + e -1 Aw = F, 

d t v - e _1 A b = G, ( 5 - 2 ) 

„ (b,v) t =o = (b 0 , v 0 ), 

where A = \/—A. Then, for any s£R and T € (0,oo] ; the following estimate holds: 


II (Mil 


Br A<7 p 2 r ) 


MIU + MMI(M)|| s+jv(4 _i )+ j 7 _ 1 , (5.3) 

24 p ' ‘ > 


with 


P> 2, - < min(l, 7 (p)), (r,p, N) / (2, oo, 3), 
r 

P> 2, ^<min(l, 7 (p)), (f,p,N) ^ (2,oo,3), 
w/iere 7 (g) := (A - 1)(± - 1), J + A = 1, and ± ^ = 1. 

Let us make the following change of functions: 

c(t,x) := eb e (e 2 t, ex), v(t,x) := eu e (e 2 t, ex), h (t,x) := d £ {e 2 t,ex). 
Then ( 6 <E ,u £ ,d e ) solves ( 11 . 111 ) if and only if (c, v, h) solves 
dtc + divv = —div(cv), 

„ /jAv + (ja + A)Vdivv P'(l + c) 

^v + v- Vv - - -M- - - 1 -A---Vc 


1 + c 

= ——div(Vh 0 Vh - -|Vh| 2 I), 
1 + c 2 

dth + v • Vh = 6>(Ah + | Vh| 2 h), 

(c, v,h)| f=0 = (co, v 0 ,h 0 ), 
with (c 0 ,v 0 ,h 0 ) := (eb £ 0 (ex), eu e 0 (ex), d £ 0 (ex)). 


1 + c 


(5.4) 
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According to Theorem ll.il there exist two positive constants rj = rj(N, A, p, £, 6, P ) and M = 

N_ 

M(N. A, p, £, 6, P ) such that (15.41) has a solution (c, v, h) in as soon as 

INI -4,oo + ll v o|| . 4 -i + ||ho — d|| n <r]. 


B, 


nW 

-° 2,1 


^ 2,1 


In addition, we have the estimate 


(c,v,h - d)|| 4 < M ||c 0 || -4 oo + ||v 0 || . + ||h 0 - d|| . n 

<b 7 v bJ' “ 


■ d 2,1 


n~2" 

■ d 2,1 


According to the scaling properties of Besov space in Proposition 12.11 it is easy to verify that 


I coll n , oo + ll v o|| - 4 -i + l|ho — d|| n « ||6g|| n <00 + 11 uq 11 n_ 1 + ||do - d|| n , 

B v ^2,1 -®2.1 ^2,1 ^2,1 


and 


(c,v,h)|| n 
<87 


, u, d) 




which combined with Proposition 15.11 conclude the part (i) of Theorem 11.21 
Now we turn to the proof of part (ii) of Theorem 11.21 

Case 1 : N > 4. Let us first focus on the convergence of ( b £ , Qu £ ). Applying Q := VA _1 div 
to the second equation of (II. lip , we conclude that ( 6 £ , Qu £ ) satisfies 
divQu £ 


d t b e + 


= — div(6 £ u £ 


d t Qu e — uAQu e + 
£ 


V6' 


= Q( — u e • Vu £ — 
e V 1 + eb e 


1, 


eb e V6 £ 

Au e + K(eb e )— 


, div(Vd e © Vd £ — -|Vd £ | 2 I)). 

1 + eb e K 2 1 V 

Setting l e := A _1 divQu e , then system (15.51) becomes 

* dtb e + e _1 AZ £ = F e 


a l e - e" 1 A6 £ = G 


(5.5) 


(5.6) 


with 


F e := —div(6 e u e ), 

/ fh e X7h e 

G e := uA/ £ - A _1 div( u £ • Vu £ + ,^ u £ - K(eb e ) - 


1 + eb e 

+ _J_div(Vd £ 0 Vd e - i|Vd e | 2 I 

Remark that Qu £ = —VA _1 Z £ so that estimating Qu £ or l e is equivalent (up to an irrelevant 
constant). Taking p = 2, r = oo, s = ^ — 1 and r = 2 in Proposition 15.21 yields 

||( 6 e , Qu e )||_ s i <Ce^m,Qnl)\\ * , + Ce^||(F £ ,G £ )|| , 




r"2 "~ 

-° 2,1 


'^(SaTi *) 


for all p > pn := .P . Note that 


||div(6 £ u £ )|| jv_ 1 < ||6 £ || jv ||u e || 

x ) lWP "lWi) 


* <C?, 
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|u £ • Vu £ 


< 11 u e II 2 „ < lln £ 


u 


V , U 


l ) \) ) mB £; ) 


< r tl/ 

+i A: l -'o > 


eb e 


l + e¥ 


Mu e 


n _i < e||6 £ || jv Mu e || i v_j 

Li(^ ) ^Wi ) 


< 


u~ l W' 


N „ U 
jT' 


) Li(S 2 y ) 


A + v < Q", 


Moreover, 


AT(eA 


W 


<ni 2 , 4 <cr 

) lhb^) 


—-—div(Vd e 0 Vd £ - -|Vd £ | 2 L 
1 + eb e K T 1 ' 


i 


1 + e& £ 


-div 


LHBZ ) 


iv(V(d £ - d) 0 V(d £ - d) - V(d £ - d)| 2 I 

< c; 


N 

^Wi ) 


<(l + e||6 e |L ,^jK-d|| 2 _..* +1 <<T 


2 

^aT*) 


' L°° (B^i) 

Collecting all the estimates above, we conclude that 


||<fc e ,Qu«)|| » . <cfd, p> 2( ” 1} . 

£Wi 7 ) N — 3 


(5.7) 


On the other hand, define A := Id — Q , then the rest part (’Pu € ,d e ) of system (11.1111 satisfies 


d f Pu £ - uAVu e = V( - u £ • Vu e -—iu f -div(Vd £ 0 Vd e )), 

1 + eo £ 1 + eb e 


(5.8) 


k d 4 d £ + u £ • Vd £ = 0{ Ad £ + |Vd £ | 2 d £ ). 

Letting w £ = Vu e — u and d £ = d £ — d, then it follows from (11.1011 and (15.81) that (w £ , d £ ) solves 


<9(W e — /iAw e = — V(u e ■ Vu c — u • Vu) — V 


eb e 


1 + eb € 


Mu* 


-V 


e 


1 + e6 £ 


div(Vd £ © Vd £ ) - £div(Vd © Vd) , 


k d t d e - 6Ad € = -u £ • Vd £ - w £ • Vd - Qu £ • Vd + d\ Vd| 2 d £ + 0\ Vd £ + Vd| | Vd £ |d £ . 


(5.9) 
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We infer from the estimates for heat equation (see ID El) that 


W e (f)|| N_ 3 +//||w e || N+ 1 


B, 


P,1 




< ||Wq|| ,jv _3 + \\V(u e ■ Vu e -u- Vu)|| 


B. 


P,1 


2 ) 


3 + 


n-^-Au- 

1 + eb e 


p‘(4- ! ) (5 - 10) 


+ 


iP 


-div(Vd e 0 Vd e ) - £div(Vd 0 Vd 


-1 + eb e 

By virtue of Corollary 12. II and interpolation, there holds 


JV_ 3 . 

LHBI i *) 


7? I ^ 


1 + e&' 


-Au ( 


' _ 3 4 


eY 


1 + e& £ 




.V _ 3 ^ , 

lHbZ 2 ) 


< ||e 6 e || .4 ||iu £ 


L3(B 2 -T *) 


^ i||6e|l ^|-*, l|u ' ll u ( 4C)~ (cr)2£i ’ 


where we have used that 


and 


_ JV i JV JV 

b G L l (B e l ’) n L°°(B e l ’ ) ^ 6 G L m {B e V ), for 1 < m < oo, 


. . jV_|_i . JV_-| I _2_ 

u G L 00 ^ ) n L\Bfc ) => U G L m (B 2 2 1 m ), for 1 < m < oo, 


(5.11) 


ll&l _^,4 ~ ll&l .f-1 + (eu)2||6 £ || . N => \\b c \\ ,N < e 2||6 C || _jv 
n? 2 ’ B 2 2 B 2 B 2 B 2 

D 2,1 • d 2,1 ■ d 2,1 


1 B 2 ’ 


Noting that 


u e • Vu e — u • Vu = u e • VQu e + u e • Vw e + Qu e • Vu + w e • Vu, 


which, together with Corollary 12.11 (15.71) . interpolation and Young inequality yields, for any 
8 > 0 

llT^u 6 • Vu e — u • Vu)|| 3 

L 1 (B p 2 ) 


< C'||u e || jv ||VQu e || jv _3 + C1|Vu|| jv ||Qu e || jv_ i 
L 2 (B^) B 2 (B p P i 2 ) L2 ( B Z ) L^B^ 2 ) 


+ C [ (||u e || jv || Vw c || jv_3 + llVull jv|| w e || N_i^\dt' 

Jo V B P i ^ B p P tl ^ 

C 5 f (||u c || 2 jv + ||u|| 2 jv ) | 

Jo v bt B a y/ 


< C'(C'o^) 2 e2 + (i||w e || jv . i + 


W JV _ 3 


B. 


i> P 2 


dtt'. 


P,1 


(5.12) 
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Moreover, 


7>(_l_div(Vd £ 0 Vd £ ) - £div(Vd © Vd)) 

£div(Vd £ © Vd e - Vd © Vd) 
eb 


N__ 3 

LHBI i 7 ) 


< 


N_ _ 3 

LHB& 2 ) 


+ 


1 + eb e 


£div(Vd £ © Vd £ ) 


N__3 

L 1 {Bl a *) 


(5.13) 


■ =h + h- 

For the estimate of I\, we have 


/i<l|Vd £ + Vd|| , ||Vd £ || V _1 

L\B^) L 2 (Bf tl ^) 


< lid 


N ,, + ||d|| JV,, 

LW, ) lWi ) 


I jV , 1 . 
L 2 (B P P i 5 ) 


For I 2 , it follows that 

/ 2 < e" ue 

4 r^i 


N _ 1 


div(V(d £ - d) © V(d £ - d)) 




(5.14) 


< 62 


N _! 


€ 1 / 


N 


“ ||d £ — d || 2 v 


+1 


l 2 {b; a ) 


(5.15) 


< 6 2 


r n || d e — d || 2 , +1 <(C?fel 

L°°{Bj ) L 2 (B 2 V ) 


Thus, substituting (|5.11|) - (|5.15|) into (15.101) . and choosing 5 in (15.121) sufficiently small, it is not 
difficult to obtain 


|w £ || jv_ 3 + -d|w £ 

* v 77 A 

'p, 1 




I JV + 1 
LHBl 1 


< IIwq|| .a 3 + Cq v e 2 + 

B P P ,i 


n e ||2 1 |||| 112 

u || . JV l || U || . N 


W £ || N_3dt' 

B p , 2 

p, 1 


(5.16) 


+ lid 


L 2 (B 2 3 +1 ) ' " "L*{B^y" "L> { B*P) 

In order to close the estimates of w £ , we now aim to bound the term ||d £ || n ,1 of ( 15.161 ). 

W/l *) 

For this, we take advantage of the estimates of heat equation (15.9D o and obtain that 


+ ||d| 


+1 


| d c 11 jv_i + 0||d £ || jv,3 

L-iB^ *) 


LHB* 




+ ||u £ • Vd £ 


B. 


p,i 


jv_i T ||w c • Vd|| n 1 T ||Qii £ • Vd|| n 1 ~\ 7 \ 

LHBfl W p P 1 2 ) W/i 3 ) (5 ' 17 ^ 


2 J£| 


+ |||Vd|^d 


LHB^ © 


_i + |||Vd £ + Vd||Vd £ |d £ 


I jv _ 1 • 
Ll (s p i 7 ) 
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Next, we estimate the terms of the right hand of (15.171) one by one. 


|u e -Vd e || jv_i < 11 u e 11 n II Vd e II n_ i 

L 2 (B£ l) 7 ) 


i 2 ) 


u 


AT 

b 2 ^) 


|d e || ^ 


||d e ll 2 

A_1 l| u II A , 3 ' 

B»(B* 2 ) 5 ) 


(5.18) 


|w e ■ Vd|| ,a_i <||w e || N-i ||Vd|| n 

- B 2 (B p p x *) 


LWi 2 ) 


< IIw e ll 2 


1 

r e \\ 2 

JV _ 3 || W || JV 


1 ll d ll 


+ 1 ^ * 




(5.19) 


||Qu e -Vd|| n_ i < || Qu e || n_i || Vd| N < || Qu e || a_i ||d|| jv + 1 . /c on) 

2 ) W/rV ^(^i) i 2 (s P p r 2 ) i a (B a Ti ) (5,20) 


|||Vd| 2 d e || n_ i < || | Vd| 2 1| a ||d e || jv_i < ||d|| 2 ^, ||d f || jsr_i . 

LHBJ i 5 ) B°°(B P T 2 ) B 2 (B 2 a ) L °°(B p P ,i 2 ) 


(5.21) 


|||Vd e + Vd||Vd e |d e 


| jV _ 1 

^(B p P ! *) 


< 


|||Vd e ||Vd e |d e || a i + |||Vd||Vd e |d e || a i 

B 1 ^ 7 ) Bi(B p 7 ) 


< 


|||Vd e |d e || w. ||Vd e || a i +|||Vd|d e || a ||Vd e || a i 

L 2 (B 2 ^) L 2 (B/i 2 ) B 2 (B.£) ^(B/i *) 


(5.22) 


< 

r-^i 


l|Vd e || 


|d e 


B 2 (B 2 ^) B»(B 2 T) B 2 (B P 5 ) 


A+i + IIVd|| a ||d ( 


JV, 1 

B 2 (B 2 3 () i»(B 2 T) L 2 (B p y 5 ) 


< lid 


'b 2 (b£ + 1 ) ' " "^(B^)/" '^“(Bjr 5 )" "^(B^ 5 ) 

It follows from the above estimates (|5.17I) - (|5.22D that 


+ ||d|| 


r+l, 


d £ 


a 11 d 
R T 


e|| 2 


l ll d ' 


l 

eII 2 


A i + ^l|d e II _ . A + 3 


B°°(B p P i 11 ) 

< K|| a _i + ||u e 
°p,i 


2 ) 




1 

iei| 2 


1 

IH e II 2 ||r] e .. 

, rt .__ ,| u || JV_1 11^ II JV i 3 

L ( B 2 i) L°°(B p p 1 7 ) L 1 ^! 2 ) 


+ IIw e || 2 JV_ 3 IIw e || 2 JV + 1 ||d|| ^ 

5 ) B 1 (B P 2 ) B 2 (B 2:i ) 


JV i-i 

"2” ' 1 n 


(5.23) 


+ IIQu'H _ A 1 ||d| 


L 2 (B p P i 2 )" "B 2 (B^i +1 ) " " L 2 { b^ 1 +L )" "L-o(B p P A ^) 


+ ||d|| 


+ l x 


+ lld 


r + l 


+ ||d| 


+i 


A ||d £ ||* 

TTi 


L 2 (B 2 y ) B 2 (B 2 y )/ I»(BJ) L oo ( fi7-5) 


1 l|d 


1 

eit 2 
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Now combining (15.161) with (I5.23[) and using Young inequality, we get 

ii *ii Mm 

W 1V_ 3 + 7T w iv.i + d n_i + — d N, 3 

L°°(B p i 2 ) 2 L i { b p 2 ) 2 ) 2 L i ( b p 2 ) 


< 


w 0 II . N_ 3 


B, 


F _ 2 

p,i 


+ || dp || n_i 
B p , 2 

P,1 


+ C^e 2 


+ / (K(0ll 2 ^+llu(0ll 2 - +lld(Oll 2 - +1 )l|w e (Oll.$ 3^' 


rT 

^* 2,1 


rT 

^* 2,1 


RT 

n 2,l 


B p 2 

P ,1 


+ / (ii<mn 2 , + i + iid(on 2 f+1 + iiu e (oii 2 f jiid^oiL#-*^. 

JO “2,1 “2,1 “2,1 


B p , 2 

p,i 


Thus, Gronwall’s inequality guarantees that 


w 


M, 


N__3 + — ||w L II JV.l + 11 d e 


L oc( B P 5 ) 2 


JV_1 + —||d e || N,3 

L °° {B p t 2 ) 2 BpB^ 2 ) 


< llwgll . iV 3 + lid, 


'0 

B p , 

p, l 


0 . A_ 

B p , 

p,i 


+ C^e 2 . 


Case 2: N = 3. Applying Proposition 15.21 to (|5.5I) with p = 2, r = oo, s = ^ and 
similar to dS3D, we obtain 

||(6 e , Qu £ )|| _ 2 p 2_1 < H(6 e ,Qu e )||_^, 2_1 <C^e2 _ r, p>2. 

L^(B p ^) L¥^(B p ^) 

Use the following interpolation for 2 < q < +oo, 


L 2 (^M + ;S^“ i 9)/(2 ‘ ?+4) ) = 


T 1 (R + ;5| 1 ); LT^(R+;S„ ? 2 ' 


2,1 / ’ ^ ’ q,l 

Make the change of parameter p = . Due to (|5.26[) , there holds 


2 

q-\- 2 


||Qu e || .4 1 < Cq U 6 2 p,2<p<+oo. 

LHB p 1 ^) 

In the following, we want to prove that 


I w e II 4 _ 3 + II w e II 4,i + ||d e || 4 _ i + ||d e || 

i°°( K ! 2 ) L'{B pTS ) l~(b p ^) li(b p ;*) 


4,3 

j- + "' 


< 


Wnll 4_3 + II dn II 4_1 TCn^e 2 P. 
U rP 2 U RP 7 ' U 


To this end, similar to the estimates of (15.1 1 D - ()5.15f) . we have 


pi ALi„. 

1 + e& £ 


-i+i Pu e 


4_3 6 t/ ii i _i_ j 

LH&l i 2 ) i”(B 2>1 5 ’)" "^(^V) 


eM- 2 -p||&T - , («/„ .. 

L°°(B^ 1 ) i) 


1,1 

3 ) 2 P||U £ || .5 


BUB 2 ) 


< 


e 2 p 


" U 


B-(B 2 ’ )" "BUS?!) 


< (C ^ 2 


T 2 p, 


(5.24) 


(5.25) 


- JiL 

p—2 ’ 


(5.26) 


(5.27) 


(5.28) 


(5.29) 
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\\V{u e - ViT-u- Vu)|| 


4 _ 3 

LHBF i 2 ) 


<C||u £ || N ||VQu e || 4_3 +C||Vu|| * ||Qu e || 4_1 

L2 (B^) L 2 (B^ x 2 ) L 2 (B^ ) L 2 (B^ 2 ) 


+ C L ( l|ue| U#J |Vwe|l B^r® + l|Vu|| siJ |w£|l Rl-0^ / 


B, 


p, i 


n 2,l 


B, 


p, 1 


(5.30) 


< C{C e 0 u ) 2 e^~p +5||w e 


Moreover, 


L 1 ^! 7 ) 


+ Cs 


u e || 2 3 + Hull 2 3 ) ||w e || 43 dt' 


B 1 - 

° 2,1 


£> 2,1 


>P 2 
J P,1 


P (rr^ div ( vd£ 0 vd ^ ~ ^ div ( vd ° vd )) 

£div(Vd e © Vd e - Vd 0 Vd) 
eb 


4 _ 3 

LHBh 


< 

r~KJ 


4 _ 3 

L 1 ^! 


+ 


l + e&‘ 

For the estimate of H\. we have 


-£div(Vd e © Vd e 


4 _ 3 := H i + H'2 ■ 
p,i ' 


L 1 ^, 2 ) 


(5.31) 


^l<l|Vd £ + Vd|| 3 ||Vd e || 4—1 

i2 (^l) i 2 (B p P ! 2 ) 


< ( ||d e || . 5 + ||d|| . 5 ||d e || 4,1 . 

L 2 ^) L 2 ^)/ L 2 ^ 2 ) 


(5.32) 


For i^ 2 i d follows that 
H 2 <e 


L°°(B 2 , /) 


div(V(d e - d) © V(d e - d)) 




e(eiy) 2 r||6 e || 2 p i (ev 
L°°{Bh) 


.3 )^ + P||d e -d|| 2 3 

L^iBlJ L "^ +1 ' 


i 2 (S P i ) 


<^-V|| 3 ||d f -d|| 2 . <(C^)M-V 

) L 2 (S|i) 

Therefore, in view of the estimates of heat equation, similar to (|5.16p . we have 


(5.33) 


|w e || 4_3 + — all W £ II 4,1 

L™{Bl i 5 ) 4 lHbp^) 


< 11 Wq 11 . 4-3 + C'q'T 2 P+ U e || 2 3 + 1 1 U11 2 3 ) || W e || 4_3 dt' 


+ ( IId' 



^i 




B l i 


5 + lldll 5 ) ||d f II 4, 1 

l 2 ^) £ 2 (B|iV l 2 !^ 2 ) 


(5.34) 
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In order to close the estimates of w £ , we must bound the term ||d £ || 4,1 of (15.341) . 

- ;p + 5 -1 


purpose, similar to (15.171) . it follows that 


L*(Bl 1 2 ) 


1 £ || 4 _ 1 + #||d £ || 4,3 

2 ) LHBl +*) 

< lldnll 4 i+||u £ -Vd £ || 4 1 + ||w £ • VdII 4 1 + ||Qu £ -Vd|| 4 1 


+ llivd| 2 d e n ,4_i + mvd- + Vd||va e |d e n . 4 - 1 . 

LHBl[, *) LHB^ 7 ) 


Next, we estimate the terms of the right hand of (15.351) as follows. 


|u £ -Vd £ || 4 _i < 11 u e 11 3 II Vd e II 4 _i 

, 7 ) L 2 (Bl 1) L 2 (B P 1 7 ) 


1 

jen 2 


< IIu e || .3 ||d £ || 2 4 _! ||d £ M , 

i2 ( B 2 2 l) L°°(B P 7 ) L!(B P J 2 ) 


|w £ • Vd|| 4 _ 1 < ||w e || 4 _ 1 ||Vd|| 3 

LHBi:*)" LHBl 1 ) 


< II W £ || 2 4_ 3 ||w e || 2 4 j ||d|| 5 , 

L°° (Bf 1 7 ) LHB p 7 ) l HB 2 - :1 ) 


||Qu e • Vd|| , 4 _i < ||Qu e || IIVd|| 3 < ||Qu e || 4 1 ||d|| 5 

- L 2 (B P X 7 ) LHB 7 ,) L 2 (B P ! 2 ) LHBlH 


'LHBl 1 2 ) 


ll|Vd| 2 d e || , 4 i < |||Vdp|| 3 ||d e || 4 i <||d 


LHBl, 2 ) 


LHBl,)" /.v/h; ; 2 ) "LHBl,)" "L°°(Bl, 7 ) 


• TT >7 . 5 


|||Vd £ + Vd||Vd £ |d £ 


L 1 ^! *) 


< 


III Vd £ | | Vd £ |d £ || 4 _i +|||Vd||Vd £ |d £ || 4 _i 


< 

r^j 


|| | Vd £ |d £ || 3 ||Vd £ || 4 _i + II | Vd|d £ || 3 ||Vd £ || 4 1 

LHBl,) LHBl, 7 ) LHBl,) LHBl, 7 ) 


< 


|| Vd e || 3 ||d e || 3 ||d e || 4 ,i + || Vd|| 3 ||d e || 3 ||d £ || 4,1 

LHBl,) L°°{Bl,) LHBlf 7 ) LHBl,) L~(B|i) L 2 (B„ P 


1 

ie || 2 
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It follows from the above estimates (I5.35l) - (|5.40p that 

||d e || 4 _ 1 +0||d £ || 4,3 

L°°(B P I 7 ) 


< lldoll 4_1 + ||u e 


A o II . 
B, 


P, 1 


L 2 (B£ i) 


|d e || ^ 


1 

ie II 2 


4 _ 1 || u || 4,3 

L°°{B?B) LHBZ i 5 ) 


+ w 


l l 

£ II 2 llw e ll 2 

4 _ 3 ll W II 4 


1 ll d ll 




+ ||Qu e || 4 _i ||d|| 5 + ||d|| 2 5 ||d' 

L 2 (Bp i 2 ) i 2 (B 2 2 i) 


L 2 (B2a) "L°°(B?J) 


+ lid 


£ 2 (4 2 i) 


5 + ||d|| 5 ) lid 6 !! 3 ||d e || 2 

L 2 (Bl i)/ 


4 1 

““■2 


| d e II 2 


L°°(B? 2 ) L 1 {Bp j 2 ) 


Now combining (15.3411 with (15.41 [) and using Young inequality, we get 
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Gronwall’s inequality then yields that 


0 u - eu 

\VT\\ 4_3 +— || w" || 4,1 + d 4_1 4- d 4,3 

L°°{Bl i 2 ) 2 2 ) ^°°(B p P i 2 ) 2 Ll (4 P i ’) 


U I, ,, 

4_3 H-W 1 


< 


Wnll 4 — 3 4“ 11 dr, 11 4_1 ) 4“CrjY 2 P. 


B p , 1 

p, 1 


B p , 2 

p, 1 


Case 3: N = 2. Applying Proposition 15.21 to (15.51) with p = 2, r = oo, s = 0 and r 
yields 
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Make the change of parameter p = 6 ^q ■ Thanks to estimate (|5.44l) , we conclude that 
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Next, we are going to prove the convergence of w £ and d £ . For this purpose, it follows from 
Corollary 12.11 that 
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Attention is now focused on bounding K\ and K 2 . For the estimate of K \, we have 
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Therefore, in view of the estimates of heat equation, similar to (I5.16p . we have 


I W e II 5 _ 5 + -u||w e || 5 ,3 

3 ) 4 LHB ^ 3 ) 


<IKII 5-1+^-^ + 

B P,1 


J u ‘ll| il + ll“ll|. il ) | |w‘IL^-|<B' 


(5.51) 


+ lid' 


L 2 (Bl ,) + H d llL 2 (B| :l) 


5 ,3 

L 2 {B% 3 ) 


In order to bound ||d e 


L 2 (B?i 3 ) 


5 _|_ 3 , similar to (15.171) . we have 


5 i +@||d c || 5 _l7 


< lldnll _5__i T 11u e • Vd e 


b pi 

+ |||Vd| 2 d e || 5 _i +|||Vd e + Vd||Vd e |d 


5 i + ||w e • Vd|| 5 i +||QiT-Vd|| 5 i ft mi 

W P T*) lhb%~*) 


LHB^i 3 ) 


L'(B% V 


Next, we estimate the terms of the right hand of (15.521) by Corollary 12.11 and Lemma 12.21 as 
follows. 
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It follows from the above estimates (I5.52I) - (I5.57I) that 
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Thus the proof of Theorem 11.21 is completed. 


□ 


References 

[1] Bahouri, H., Chemin, J.-Y., Danchin, R.: Fourier analysis and nonlinear partial differential equations, volume 
343 of Grundlehren der Mathematischen Wissenschaften. Springer, Heidelberg, (2011) 

[2] Bony, J.-M.: Calcul symbolique et propagation des singularity pour les equations aux derivees partielles non 
lineaires. Ann. Sci. Ecole Norm. Sup. (4) 14 (2), 209-246 (1981) 

[3] Chemin, J.-Y., Gallagher, I., Iftimie, D., Ball, J., Welsh, D.: Perfect incompressible fluids. Clarendon Press 
Oxford, (1998) 

[4] Chemin, J.-Y., Lerner, N.: Flot de champs de vecteurs non lipschitziens et equations de Navier-Stokes. J. 
Differential Equations 121 ( 2 ), 314-328 (1995) 

[5] Danchin, R.: Global existence in critical spaces for compressible Navier-Stokes equations. Invent. Math. 
141 ( 3 ), 579-614 (2000) 

[6] Danchin, R.: Zero Mach number limit in critical spaces for compressible Navier-Stokes equations. Ann. Sci. 
Ecole Norm. Sup. (4) 35 (1), 27-75 (2002) 

[7] Desjardins, B., Grenier, E.: Low Mach number limit of viscous compressible flows in the whole space. R. 
Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci. 455 (1986), 2271-2279 (1999) 

[8] Ding, S., Huang, J., Wen, H., Zi, R.: Incompressible limit of the compressible nematic liquid crystal flow. J. 
Funct. Anal. 264 (7), 1711-1756 (2013) 

[9] Ding, S., Lin, J., Wang, C., Wen, H.: Compressible hydrodynamic flow of liquid crystals in 1-D. Discrete 
Contin. Dyn. Syst. 32 (2), 539-563 (2012) 

[10] Ding, S., Wang, C., Wen, H.: Weak solution to compressible hydrodynamic flow of liquid crystals in dimension 
one. Discrete Contin. Dyn. Syst. Ser. B 15 (2), 357-371 (2011) 

[11] Ericksen, J. L.: Conservation laws for liquid crystals. Trans. Soc. Rheology 5, 23-34 (1961) 

[12] Ericksen, J. L.: Hydrostatic theory of liquid crystals. Arch. Rational Mech. Anal. 9, 371-378 (1962) 

[13] Fang, D., Zi, R.: Incompressible limit of Oldroyd-B fluids in the whole space. J. Differential Equations 256 
(7), 2559-2602 (2014) 

[14] Ginibre, J., Velo, G.: Generalized Strichartz inequalities for the wave equation. J. Funct. Anal. 133 (1), 
50-68 (1995) 

[15] Hao, Y., Liu, X.: Incompressible limit of a compressible liquid crystals system. Acta Math. Sci. Ser. B Engl. 
Ed. 33 (3), 781 -796 (2013) 

[16] Hong, M.-C.: Global existence of solutions of the simplified Ericksen-Leslie system in dimension two. Calc. 
Var. Partial Differential Equations 40 (1-2), 15-36 (2011) 

[17] Hu, X., Wu, H.: Global solution to the three-dimensional compressible flow of liquid crystals. SIAM J. Math. 
Anal. 45 (5), 2678-2699 (2013) 

[18] Huang, T., Wang, C., Wen, H.: Blow up criterion for compressible nematic liquid crystal flows in dimension 
three. Arch. Ration. Mech. Anal. 204 (1), 285-311 (2012) 

[19] Huang, T., Wang, C., Wen, H.: Strong solutions of the compressible nematic liquid crystal flow. J. Differential 
Equations 252 (3), 2222-2265 (2012) 

[20] Jiang, F., Jiang, S., Wang, D.: On multi-dimensional compressible flows of nematic liquid crystals with large 
initial energy in a bounded domain. J. Funct. Anal. 265 (12), 3369-3397 (2013) 

[21] Jiang, F., Jiang, S., Wang, D.: Global weak solutions to the equations of compressible flow of nematic liquid 
crystals in two dimensions. Arch. Ration. Mech. Anal. 214 (2), 403-451 (2014) 


36 


Q. BIE, H. CUI, Q. WANG, AND Z.-A. YAO 


[22] Jiang, S., Ju, Q., Li, F.: Incompressible limit of the compressible magnetohydrodynamic equations with 
periodic boundary conditions. Comm. Math. Phys. 297 (2), 371-400 (2010) 

[23] Keel, M., Tao, T.: Endpoint Strichartz estimates. Amer. J. Math. 120 (5), 955-980 (1998) 

[24] Leslie, F. M.: Some constitutive equations for liquid crystals. Arch. Rational Mech. Anal. 28 (4), 265-283 
(1968) 

[25] Lin, F.: Nonlinear theory of defects in nematic liquid crystals; phase transition and flow phenomena. Comm. 
Pure Appl. Math. 42 (6), 789-814 (1989) 

[26] Lin, F., Lin, J., Wang, C.: Liquid crystal flows in two dimensions. Arch. Ration. Mech. Anal. 197 (1), 
297-336 (2010) 

[27] Lin, F., Liu, C.: Nonparabolic dissipative systems modeling the flow of liquid crystals. Comm. Pure Appl. 
Math. 48 (5), 501-537 (1995) 

[28] Lin, F., Wang, C.: Global existence of weak solutions of the nematic liquid crystal flow in dimension three. 
Comm. Pure Appl. Math. (To appear) 

[29] Lin, J., Lai, B., Wang, C.: Global finite energy weak solutions to the compressible nematic liquid crystal 
flow in dimension three, http://arxiv.org/abs/1408.4149, arXiv:1408.4149 

[30] Lions, P.-L.: Mathematical topics in fluid mechanics. Vol. 2. Compressible models, volume 10 of Oxford Lec¬ 
ture Series in Mathematics and its Applications, Oxford Science Publications. The Clarendon Press, Oxford 
University Press, New York, (1998) 

[31] Lions, P.-L., Masmoudi, N.: Incompressible limit for a viscous compressible fluid. Journal de mathematiques 
pures et appliquees, 77 (6), 585-627 (1998) 

[32] Ou, Y.: Low mach number limit of viscous polytropic fluid flows. J. Differential Equations 251 (8), 2037-2065 

( 2011 ) 

[33] Runst, T., Sickel, W.: Sobolev spaces of fractional order, Nemytskij operators, and nonlinear partial differ¬ 
ential equations, volume 3. Walter de Gruyter, (1996) 

[34] Strichartz, R. S.: Restrictions of Fourier transforms to quadratic surfaces and decay of solutions of wave 
equations. Duke Math. J. 44 (3), 705-714 (1977) 

[35] Wang, D., Yu, C.: Incompressible limit for the compressible flow of liquid crystals. J. Math. Fluid Mech. 16 
(4), 771-786 (2014) 

[36] Xu, F., Hao, S., Yuan, J.: Well-posedness for the density-dependent incompressible flow of liquid crystals. 
Math. Meth. Appl. Sci. 38 (13), 2680-2702 (2015) 

College of Science, China Three Gorges University, Yichang 443002, PR China 

E-mail address: qybie@126.com 

School of Mathematical Sciences, Huaqiao University, Quanzhou 362021, PR China 

E-mail address: cuihaibo2000@163.com 

School of Mathematics and Computational Science, Sun Yat-Sen University, Guangzhou 510275, 

PR China 

E-mail address: mcswqr@mail.sysu.edu.cn 

School of Mathematics and Computational Science, Sun Yat-Sen University, Guangzhou 510275, 

PR China 

E-mail address: mcsyao@mail.sysu.edu.cn 


